Derivation of hydrodynamics for the gapless mode in the BEC-BCS crossover from 

the exact one-loop effective action 
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We show that many hydrodynamical properties of the BEC/BCS crossover in the presence of a 
Feshbach resonance at T — can be derived easily from the derivative expansion of the (exact) fully 
renormalized one-loop effective action. In particular, we calculate the velocity of sound through- 
out the BCS and BEC regimes and derive the generalized superfluid continuity equations for the 
composite two -fluid system. 
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As the strength of fermionic pairing increases in cold 
alkali atoms there is a continuous evolution from the 
BCS-like behavior of Cooper pairs to Bose-Einstein Con- 
densation (BEC) of molecules. This crossover phe- 
nomenon becomes controllable experimentally [l I when 
the atoms interact through a Feshbach resonance 1, d, 0, 
01- In both regimes, condensation is a consequence 
of global U(l) symmetry breaking with a concomitant 
gapless Goldstone mode (the phonon). 

Although the s-wave scattering length diverges at the 
crossover, the speed of sound changes smoothly, as do the 
hydrodynamical properties of the condensate. The main 
result of this letter is to show that, at T — 0, all of these 
can be determined straightforwardly from the derivative 
expansion of the exact, fully renormalized one-loop effec- 
tive action, in the spirit of [7| , sidestepping the more com- 
plicated derivations based on a full multi-channel analysis 
(e.g. see 0|). 

Consider a condensate comprising a mixture of 
fermionic atoms and molecular bosons. The fermions 
ip a (x), with spin a = (t,|), undergo self-interaction 
through an s-wavc BCS-type term. In addition, two 
fermions can be bound into a molecular boson <j)(x) 
through a Feshbach resonance effect. To exemplify our 
method, we take the Lagrangian density to be 0, Q 
{U > 0, g fixed) 



+ u r^x) $i{x) ^o) ^(x) 



V 2 

i d t + - h 

Ira 



- g [4>*(x) ^(x) ^(x) + c)>(x) $l(x) V T »] • (1) 



The mass of the bosons is twice that of the fermions, M — 
2m. The kinetic energy of the fermions is = k 2 /2m 



and the kinetic energy of the bosons is k 2 /2M + v, where 
v is the threshold energy of the Feshbach resonance. To 
demonstrate the method we restrict ourselves to a nar- 
row resonance approximation 0, 0| where quantum loo 
effects from the molecular boson can be safely ignored 
although this is often an idealization. The effective four- 
fermion interaction is determined by the value of v, which 
is controllable by an external magnetic field. On tuning 
the field, the condensate can be varied in form from BCS 
Cooper pairs to BEC molecules as the scattering length 
changes sign. 

On introducing the auxiliary field A(x) = 
Uipi(x)ip-\(x), a Hubbard-Stratonovich transforma- 
tion leads to an effective Lagrangian density, to which 
the fermionic contribution is ^ (x)G~ 1, ^(x), where ^(x) 
is the Nambu spinor, ^(x) — (i/>| , i/'j. ) , and G _1 is the 
inverse Nambu Green function 



G = 



id t - e 
A* (*) 



A(aO 
id t + £ 



(2) 



The combined condensate of the 



with e = — J u. 

theory is A(x), given in terms of the bifermion and 
molecular condensates A and <fi as A(x) = A(x) — 
g 4>{x). The gapless mode of the theory is encoded in 
the phases of A(x) and 4>{x) for which we write A(x) = 
|A(x)| e t0&( - x \ 4>{x) = -\4>(x)\ e l0 ^ x l The amplitude 
and phase of A(x) = |A(x)| e l9 ^ x ^ can be determined 
from those of A(x) and <j){x) by its definition above. We 
now perform a ?7(1) gauge transformation on the fermion 



field Vcr(x) 



je A (x)/2 



X<t(x). Integrating out x<r(x) 



leads to a non-local effective action S s[(j>, 4>* , A, A*]. 

The action possesses a U(l) invariance under the phase 
change 9^ — * 9^ + a and 9a — * 9\ + a (or 9^ — * 9^ + 
a). This symmetry is spontaneously broken when A and 
4> respectively acquire the non-vanishing constant values 
Aq and 00 determined by the gap equations obtained 



2 



from extremizing the effective action. In these, |Ao| 
|Ao| + I — g0o| satisfies 



1 



d 3 p 1 



(2tt) 3 2E p 



(3) 



with U cS = U + g 2 /(v - 2/i) and £ p = {e 2 p + |A | 2 ) 1/2 , 
and 



I0o] 



.9 



|Ao| 



(4) 



i/ - 2/i U cS 

Since the microscopic theory under consideration is 
Galilean invariant, any effective theory derived from it 
must respect this symmetry Q. Consider the dynam- 
ics of the phonon carried in the angular variables. To 
preserve Galilean invariance at each step, the varia- 
tions in the condensate magnitudes can be written as 
|A| = |A |+S|A|, \<f>\ = IM+Sm BE El. We 
assume that terms in <5|A|, <5|<^|, and (0 A — 9^,) 2 are of 
the same order in their defining equation, although #a 
and ()<), are large variables. We now use the fact that 

e -iv 3 8 K {x)/2Q-l el a 3 e^(x)/2 = Q-l _ ^ where 

S = (-iV 2 6» A /4m+ (V0 A )(-iV)/2m)I (5) 
+ [6-J2 + (V9 A ) 2 /8m)a 3 - 5|A| a u 

with 5\A\ defined by |A| = |A | + S\A\ [9\. G 1 is the 
free inverse Nambu Green function with the same form 
as G _1 in Eq.([2]), except that A(x) is replaced by |Ao|. 
S e s [4>, 4>* , A, A*] then permits the derivative expansion 

ScS = -iTrLniG^+iTrjr^^- /W A " 



n=l 
V72 



u 



+ J d*x<b*(x)(id t + — + 2 f i-vj^x), (6) 

the first two terms of which are no more than the ex- 
pansion of — iTrLn(C7 ). For our purposes it is suf- 
ficient to truncate the sum in n at second order, to 

( 2") ~ 

yield an action S^, s [S\(j)\,5\A\,9^,9 A } after eliminating 
A in favor of A, and 9 a in terms of 6 A , on using 
(0 A - 6$) fa (|Aq|/|A|)(0a - fy,)- Diagrammatically, this 
amounts to taking account of fermionic one-loop effects 
in the effective action. 

The action of the gapless phonon mode is that part of 

(2) 

the quadratic contribution to Sg ff in which the deriva- 
tives of 8\<j)\ and 5|A| are omitted. After straightforward 
manipulations this takes the form 



5, 



phonon 



l 



l 



In 2 

2 

(2/i 



2\^\5\<j>\ 04,-ae L 6\h\ 



u 



5\A\8\<t>\ 



(7) 



This lends itself to a very simple mechanical picture of 
a coupled 'wheel' and 'axle'. The radius of the wheel is 
|A| and that of the axle is | — g<j>\, measured from the gap 
values | A 1 and | — g (f>o\ with angles of displacement 9 A , 
6*0, respectively. There is slippage between the wheel and 
axle with an elastic restoring force f2 2 = (2g/C/)|0o||Ao|. 

The fermion number density at n — 1 is 



Po = Po +Po, 



(8) 



where p$ — Jd 3 p/(2n) 3 [l — e p /E p ] is the ex- 
plicit fermion density, and Pq = 2|0o| 2 is due to 
molecules (two fermions per molecule). The other 
coefficients are straightforwardly derived as N — 
/d 3 p/(27r) 3 (|Ao| 2 /2^ 3 ),a = / d 3 p/(27r) 3 (|A | £p /2£: 3 ), 
and / 3 = /d 3 p/(2^) 3 ( £ 2/i ? 3 ) . 

A further straightforward eigenvalue calculation gives 
the long wavelength dispersion relation for the phonon 
as lu 2 — v 2 k 2 + 0(k 4 ), where 



Po/2m 
N + A/B 



(9) 



and 
.4 

B 



a 2 (2^ - v) 



U 
1 / 2 



2a|. 



2 V U 



1 U 



- P (2/x - v) 



UcS 

U ' 



(10) 



Note that v is independent of the slippage strength f2 ^ 0. 

In fact, (3 is UV-singular, as are U e ff,U and g from 
the gap equations ([3]) and (j4]). Renormalization is im- 
plemented by defining the renormalized coupling U e s in 
terms of the s-wave scattering length a s as 



N 
k F a s 



1 



1 



d 3 p 1 



(2tt)3 2e F 



f A d 3 p 


1 1 " 


1 (2^) 3 


2 1p ^ 



(11) 



where kp is the Fermi momentum and A is a UV cutoff. 
UV-finite renormalized couplings U and g are defined 
similarly in the limit A — > oo ; 

1 _ 1 f A d 3 p 1 1 _ 1 U f A d 3 p 1 
U~U~J (2^2^' ~§~~g~~g~J V 

(12) 

In turn, we now define a renormalized threshold energy 
v through 



t/eS = U + 



9 



v — 2/i 



(13) 



The outcome of this renormalization is that S'phonon of ([7|) 
and the gap equation for the condensate ^ are rendered 
UV finite term by term by replacing the unrenormalized 



3 




-3-2-10 1 2 3 4 
v/2e F 

FIG. 1: The behavior of the sound velocity v throughout 
BEC and BCS regimes as a function of the threshold energy 
v. We choose U = 7.54eir/k% and g = 4.62 ep/k 3 J 2 as an ex- 
ample [T2]. The dotted line is obtained with the approximate 
solution of the sound velocity, Eq. (|15p for the BEC regime and 
Eq. (|14p for the BCS regime. The lower inset shows how the 
scattering length a s varies from BEC to BCS as the threshold 
energy increases, while the upper inset reveals the evolution 
of the fermion density p$ (solid line) and the molecule density 
po (dotted line). 



U, U e g, g, (3, v by their renormalized counterparts. Hence- 
forth we drop the overbars for simplicity, and understand 
all quantities in j9| as renormalized. Necessarily, the 
renormalization prescription we propose above makes not 
only the gap equations free of UV-divergence, but also the 
sound velocity. Thus, combining the renormalized equa- 
tions ([4| and (fTTj) with the number equation J5]) allows us 
to study the behavior of the condensate numerically when 
the threshold energy varies so that the scattering length 
goes from a s — ► + , the deep BEC regime, to a s — > 0~, 
that of deep BCS, through a BCS/BEC crossover. This 
is shown in the inset in Fig. ([1]) 0, [lj| • Solving first for 
the chemical potential /i as a function of the threshold 
energy v, the sound velocity ([9]) can be computed, as 
depicted in the main Fig.([T]), in agreement with (l3j |. 

To understand the numerical behavior shown in Fig. ([T]) 
we evaluate v 2 analytically in both deep BCS and BEC 
regimes. In the deep BCS regime, there is consider- 
able simplification as all momentum integrals are dom- 
inated by the k modes near the Fermi surface. In this 
limit, 1 0o | ~ 0, where few molecules are present, and 
a w 0, the ignorable amplitude-phase coupling due to 
the particle-hole symmetry near the Fermi surface. We 
obtain A/B <C N, leading to 



Po 
2mN' 



(14) 



With the total fermion density po = k F /3ir 2 and N w 
No = mkp/2-K 2 , the fermion density of states at the 
Fermi surface, the sound velocity is given by v 2 ~ vp/2>. 
This is the result we obtain from conventional BCS the- 
ory with no Feshbach resonance (e.g. see Q), for which 
5 = and \c/>q\ = 0. 

On the other hand, in the deep BEC limit U <C 
g 2 j\v — 2/i|, which gives rise to \v — 2p\ g 2 N /kpa s , as 
obtained from £/ e ff. To maintain this relation for small 
a s , U cannot be too large. From Eqs. (fTO]) . the sound 
velocity v 2 can be approximated as 



8m \fi\ 



(15) 



in this deep BEC regime as a result of A/B N. We 
find that \4>o\ increases as v goes increasingly negative, 
while the combined condensate | Aq| decreases to keep the 
number of fermions fixed. As a result |Ao| 2 ~ <7 2 |</>o| 2 — 
(g 2 /2)po for largely negative v (i.e. a s — > + ) when all 
fermions are in the form of molecules as seen in the inset 
in Fig.(l). Thus, the behavior of v 2 , which is found to 
approach zero at small a g , is now determined by how the 
chemical potential p increases negatively as a s deceases. 

In the central region, the behavior is smooth across the 
'unitarity limit' at \a s \ — > oo when v = 2p. In this case, 
the gap and number equations are reduced to involving 
only one parameter g in the one-channel formulation of 
Ref. [14| . If g is large, 'universal behavior' is found where 
v 2 sa 0.2vp. Nevertheless, the renormalization of the 
molecular boson is expected to contribute sizable cor- 
rections to the sound velocity obtained above in such a 
strongly coupled regime (isT ]. 

More generally, the BEC/BCS system permits a hydro- 
dynamic interpretation as a two-component superfluid. 
The crucial ingredient is the Galilean invariance of the 
derivative expansion. It is not difficult to rederive the 

(2) 

relevant angular parts of S^. s from S'phonon- We restore 
Galilean invariance by the substitutions 

• • , (ve) 2 (vef . , (w) 2 



4m 



4m 



(16) 



4m 

for both phase angles. 

(2) 

Taking the variation of S^ s with respect to the con- 
densate phase 9^ leads to a single mean-field equation of 
motion which can be rewritten in terms of the explicit 
fermion number density pf'- 

d iP F + v ■■ iF -2n 2 {e K -e 4> ) = o, (17) 



where 



PF 



P% 



N 



4m 



2a 5 A 



(18) 



and jp — pp\76^/2m. There is a similar equation for 
the fermion number density due to molecules: 

d 



df 



PB + V • ] B 



2n 2 ( 



0, 



(19) 



4 



where (to lowest order) 

p B =p° ! +4|0o|^ = 2|0| 2 , 



(20) 



and j B = pB^/O^/Am. Putting these together gives the 
continuity equation for total fermion number 



d 

w 



V-j = 0, 



(21) 



describing a coupled two-component superfluid, where 
p = pp + pb and j = j F + js- The behaviour of the 
mean explicit and molecular fermion number densities 
p F and p B as v varies is given in the inset of Fig.l. 

Our earlier definitions give fl 2 (8 A — 0$) oc |Ao||0o|(#a — 
6$), showing that the coupling between the superfluid 
components is due to the difference in the phase of the 
fermionic pairs A and the molecular field <j>. It vanishes 
in both the deep BCS and BEC regimes, when \(j> \ and 
|Aq| tend to zero respectively. Further, with Pf ^ Pb 
and 3> Us | in the deep BCS regime and pp <C ps 
and <C \ib\ in the deep BEC regime, the system 
is described by a single fluid in each case. Away from 
these extremes the situation gets more complicated, with 
the coupling strongest in the transition regime, but still 
tractable for vortices with phase coupling 9a = 0$, whose 
properties will be pursued elsewhere. 

Finally, we briefly consider Gross-Pitaevskii (GP) 
equations in the BEC regime. (For the BCS regime the 
results of @ can be simply extended to j ^ 0). The 
relationship of superfluid equations to GP equations is 
well established, in principle. To each (p,j) pair there 
is allocated a complex GP (or non-linear Shroedinger) 
field 'S. In our case the Gross-Pitaevskii fields underly- 
ing (TJI]| and (19]) are * F = y/p~F~e ie * /\/2 and * B = <j) = 
y/PB~e ie 't' I \[2. Although, in general, ^f has the phase 
of the combined condensate, but the magnitude due to 
the explicit fermion density only, it happens that, in the 
BEC regime, 



I*; 



2 _ 1 (2mlMl) 3/2 |A|2 
1 32it Id 2 1 1 ' 



(22) 



That is, now fj? a A, linking its phase to the condensate 
density. 

Furthermore, S'phonon is all that is needed to extract 
the coupling constant of two-body interactions between 
the condensate. The form of this two-body interaction is 

A A (At A) 2 = A A (|A | 4 + 4|A | 2 (J|A|) 2 + ...) (23) 

From the second term we can read off A^ directly from 
S'phonon as A A = (2/U - /3)/S| A | 2 . Again in the BEC 
regime, we find 



1 {2m\p\f' 2 
256tt \pf ' 



(24) 

Let us rewrite A A (A^A) 2 as a GP self-interaction 
AI'J'fI 4 , using the renormalization of (f2"2")) . On using 



we find a weakly repulsive interaction A = — 2ir\a A \/M, 
where a A — 2a s and M = 2m, as follows from the strong- 
coupling Bogoliubov-de Gennes equations [iH ]. 

The extension of our approach to non-zero temperature 
is straightforward, in principle, and will be considered 
elsewhere. It has yet to be seen whether, in general, the 
effect of Landau damping can be interpreted as a normal 
fluid component in addition to the cou pled superfluids of 
([2"Tj) . as happens for pure BCS theory [13] . 
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